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The experimental procedures and data analysis methods used in the determination of 
total atmospheric refraction by radiometric means are described. The results of five months' 
observation are presented in plots of the mean refraction, its standard deviation, and stand- 
ard error of estimate for specified altitude angles between 2 and 65 degrees. The a.m. values 
of refraction are significantly greater than those of the p.m. at the same altitude angle. This 
effect is attributed to the diurnal cycle. The measured total refraction exhibits a strong 
linear correlation with surface refractivity 



1. Introduction 

Development of a precise radiometric sextant has 
made possible the measurement of the effect of the 
earth's atmosphere on the direction of propagation 
of 2-cm wavelength radio waves. In particular, this 
paper considers the problem from a phenomenological 
point of view, and reports the performance and anal- 
ysis of empirical observations over a wide range of 
observed altitude angles and for the variety of 
weather conditions found in nature. Refractive 
bending effects in the vertical plane (that is, the alti- 
tude coordinate) are considered since, under the usual 
assumption of spherical stratification, azimuthal re- 
fraction does not exist. Further, ionospheric effects 
are assumed negligible at 2-cm wavelengths; there- 
fore, only tropospheric phenomena are being con- 
sidered. 

Smart [1956] has shown that, in general, atmos- 
pheric refraction depends on the index of refraction 
profile. However, if the observed altitude is re- 
stricted to relatively high angles (that is, greater 
than 20 deg), total atmospheric refraction is found to 
depend only on the surface value of refractivity and 
the observed altitude angle. Extension of the ana- 
lytical consideration to low-altitude angles becomes 
very complex, and generally ray- tracing techniques 
have been followed in this case. The usual procedure 
is to measure the refractivity profile with a radio- 
sonde instrument and then compute, from this profile, 
the total atmospheric refraction by numerical inte- 
gration methods [Shulkin, 1952]. 

In spite of the wide use of computed values of at- 
mospheric refraction, it appears that little informa- 
tion is available on precise measurements of this 
quantity at radio wavelengths. Anderson et al. 
[1960] used a beacon tracking radar to measure low- 
angle refraction through part of the atmosphere. 
Experimental studies of very low-angle (horizontal) 
refraction through the entire atmosphere have been 
performed by Aarons et al. [1958]. 

It is believed that the investigation reported 
herein represents the first attempt to precisely meas- 



ure total atmospheric microwave refraction over a 
wide range of altitude angles and with the normal 
variation of meteorological parameters of the atmos- 
phere. To the knowledge of the author, the only 
studies which closely approach the present one were 
those of Marner and Ringcen [1956], and Marner 
and Stewart [1955] at 8.7 mm. Ringoen et al, [1952] 
did perform some measurements at 1.91 cm with a 
radio sextant of significantly lower precision than 
the one used in this study. However, the refinements 
of equipment and analysis technique used here are 
believed important enough to make this study 
unique. 

The plan was to use a radiometric sextant to track 
the sun as it traversed the sky, and to record the 
observed values of the sun's altitude angle at certain 
precisely determined times throughout the day. The 
measurements were obtained at Cedar Rapids, Iowa, 
between August and December 1959. The observed 
altitude angles were limited by ground effects and 
by the maximum observed solar altitude angle to 
within 2 and 65 deg respectively. The U.S. Naval 
Observatory made available highly accurate ephem- 
eris data from which the sun's position was computed 
for any specific time for which a radio sextant obser- 
vation was made. These computed values were con- 
ceived to be true values, and any deviation from 
them was due to error-producing causes in the equip- 
ment, the atmosphere, or the center of microwave 
radiation from the celestial source. The difference 
between the observed and true altitude angle was 
the total error produced by all causes, among them 
atmospheric refraction. Fortunately, atmospheric 
refraction was the largest of all errors, so its effect 
was distinguished easily. 

Operation of a radio sextant system depends on 
microwave radiation emanating from the sun. In 
the 2-cm microwave region, the sun radiates thermal 
noise energy similar to a blackbody of temperature 
approximately 7,000 °K. The space surrounding the 
sun also radiates a very minute amount of energy 
and appears as a blackbody of temperature very 
near absolute zero. 
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A sensitive microwave receiver connected to a 
highly directional antenna senses when the sun is in 
the antenna beam by an increase in noise output 
from the receiver. If the sun were directly centered 
in the beam, the noise output would be maximum, 
and if the sun were displaced slightly from boresight, 
the output would be reduced. Assuming symmetry 
of the antenna beam and of the temperature distri- 
bution of the source, circular scanning of the sun's 
disk with the antenna beam can provide tracking 
error information. That is, if the sun is centered in 
the scan pattern, no change in d-c output level is 
obtained over the scan cycle, since the beam always 
sees the same effective source temperature. If the 
sun is not centered on the scan axis, the d-c output 
is modulated at the scan rate and correlated (phased) 
with the direction of the scan-axis to source-center 
displacement. The modulation arises since during 
part of the scan cycle the beam predominantly sees 
the cold space, and during the remaining part of the 
cycle it sees the hot sun. 

After suitable a-c amplification, the modulated 
signal is synchronously detected to provide a pair of 
d-c voltages, one in up-down, the other in right-left 
coordinates, which are proportional in magnitude 
and polarity to the magnitude and direction of the 
sun's displacement in line-of-sight coordinates. These 
signals are used to achieve automatic tracking. 

An air-supported radome was used to provide a 
temperature controlled environment for the sextant 
and to protect it from the effects of weather and 
wind. This radome was made by cementing sections 
of neoprene-coated nylon fabric into the shape of a 
truncated sphere. It was fastened at its base to a 
circular ring which was free to rotate in the hori- 
zontal plane. The ring was driven by an electric 
motor so that the radome rotated 360 deg in azi- 
muth every 20 sec. Radome rotation was necessary 
to remove refractive errors produced by localized 
differences in the radome material. Rotating the 
radome caused rapid fluctuations in apparent posi- 
tion of the source which were effectively filtered, or 
averaged by the tracking-servo time constant. 
Thus, the sextant followed the mean apparent posi- 
tion of the source. There is adequate evidence 
beyond the scope of this paper that radome rotation 
successfully removed the effects of radome structural 
variations for altitude angles greater than 15 deg. 
Below this angle the effects are somewhat uncertain, 
and this analysis will consider them. 

2. Data Acquisition and Reduction 
We have the following mathematical model: 



h=h+T,FHt(x,y, . . .)+Rn(a,b, . . .), (1) 

i = l 

where h is the true altitude angle, and h is the ob- 
served altitude angle. F hi (x,y, . . .) are functional 
relations which describe predictable (systematic 
and assignable) errors in the altitude coordinate, 
that is, those errors that are assignable to known 



causes, and are known functions of the several 
independent variables x,y, . . . and perhaps the de- 
pendent variable h. R h (a,b, . . .) is the residual error 
in the altitude coordinate, and is some unknown 
function of the infinite variety of independent 
variables a,b, . . . . That is, R h (a,b, . . .) is random 
in that we cannot predict its precise value from a 
knowledge of only a finite number of independent 
variables. 

To arrive at a model which is useful for the digital 
computations performed in this analysis (1) must be 
expressed in discrete, rather than analog form. 
Thus, for thejth observation, we have 



h j =h 0j +Ah lj -\-Ah 2j + . . . Ah kj +e j} 



(2) 



where Ah kj is the error due to the kih assignable 
cause under the specific conditions existing at the 
instant the jth observation was made. The residual 
random error at the instant of observation is e i% 
In effect then, Ah kj is a correction term due to the 
kth cause which must be added to the observed value 
of altitude angle to make it more nearly correct. 
For example, AA U might be the correction due to 
atmospheric refraction and Ah 2i might be the cor- 
rection due to a systematic error in the angular 
readout mechanism. For this analysis, the residual 
error is assumed to be normally distributed with 
mean value (/*) of zero and standard deviation, a. 

Let the first correction term, Ah u be the total 
measured atmospheric refraction, r. Then, from 
the true and observed altitude angles and the re- 
maining correction terms, the atmospheric refraction 
can be computed from \ 



Tj+ej^hj—hoj—Ah^j—. . . Ah kj . 



(3) 



As long as the conditions assumed for Bj hold, (3) 
is an unbiased estimate of the true value of measured 
refraction. Notice that refraction cannot be deter- 
mined perfectly, it can only be estimated to within 
a standard error of <r. 

The problem in making precise determinations of 
atmospheric refraction by this method reduces to 
one of evaluating the various errors introduced by 
the measuring instrument and by the sun. 

An angular error can be introduced into the 
measurement of solar altitude angle if the distribu- 
tion of intensity of microwave radiation from the 
sun should become asymmetric with respect to its 
center. Such asymmetry can result from enhanced 
radiation in localized regions of flare or sunspot 
activity. The time and location of the occurrence 
of an active region on the sun's disk was not pre- 
dictable, but the fact of its occurrence was detectable 
with the radio sextant and by other independent 
observations. Thus, intervals when dynamic solar 
activity effects were present could be identified and 
were rejected from the sample of observations used 
for the determination of refraction. Fortunately, 
the frequency and duration of solar activity at the 
observational wavelength were small and relatively 
little data were lost. 
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Ideally, it would be desirable to determine the 
instrument error corrections by independent tech- 
niques. Unfortunately no such determination of 
the readout dial-index error could be implemented 
which would have the desired accuracy. Therefore, 
it was necessary to evaluate this error by first 
eliminating all other errors from the data. The 
dial-index error is independent of altitude angle so 
evaluation of the remaining errors was required 
only over a small range of angles or, at the least, 
some specific altitude angle. 

Fortunately, total atmospheric refraction can 
be accurately determined for high-altitude angles 
by computation requiring only surface measured 
information. This permitted removal of the largest 
of the error components and revealed the remaining 
small instrument errors which were determinable 
by independent means. After the various error 
components were evaluated, the process was reversed 
and (3) was used to compute the total refraction 
for all angles of observation. 

It was found that there were three significant 
sources of tracking errors within the instrument 
itself. They are listed and discussed below. 

(1) Outer-loop bias error — a nonconstant, uncon- 
trolled, assignable error arising fro in residual biases 
in the outer-loop of the tracking-servo. 'Phis error 
was considered to be constant over various prede- 
termined time intervals and varied from interval 
to interval. Its magnitude was determined from 
the synchronous detector output recordings. Thai 
is, the output voltage recordings were examined 
to determine both the time interval over which the 
output voltage was considered to have a given 
average value, and the value of that average. The 
angular offset was calculated by dividing this value 
by the voltage-angle system analog. This analog- 
was determined periodically throughout the entire 
period of tracking activity. The angular outer-loop 
offset errors, so determined, were used to correct the 
sextant's indicated altitude angle. 

(2) Readout transmission error — a sinusoidal 
function of observed altitude angle which arises 
from slight differences in the gain and phase shift 
of the 10-kc transmission line between the altitude 
axis angle sensing synchro and its follow-up resolver. 
This error was known to be of the form A sin 
(360 A +</>). Over the time interval between adjust- 
ments of the gain and phase controls it was specified 
by determining the constants A and <j>. The phase- 
angle constant was estimated by observation of 
graphical plots. A linear correlation of the residual 
tracking errors and the mathematical function 
sin (360 /^o+0) was performed. Several trials were 
made in which different estimated values of the 
phase angle were used. The phase angle yielding 
the largest correlation coefficient was considered 
the proper one to specify the readout transmission 
error. The amplitude constant, A, was given by 
the slope of the linear regression line. The value 
of this error was computed for every tracking 
observation and the indicated altitude was corrected 
by this value. 



(3) Dial index error — ft constant value arising* 
from error in setting the readout dial with respect 
to the boresight axis. This error is considered to be 
fixed throughout the intervals between adjustment 
of the dials and/or mechanical failure or adjustment 
which would change the relationship between the 
dial and the boresight axis. This error was deter- 
mined by computing the average tracking error 
(corrected for calculated high-angle refraction and 
outer-loop error) throughout the entire interval, the 
averaging being done for all observed altitude angles 
greater than 20 deg. Again, each tracking observa- 
tion was corrected by the appropriate dial index 
value. 

Throughout the tracking period there were certain 
contaminating influences which further affected the 
measurements. These contaminations were of such 
a nature that they did not readily permit determina- 
tion of the amount of error and therefore the data 
could not be corrected with any degree of accuracy. 
Observations known to have been mad( k under the 
influence of these uncorrectable effects were rejected 
from further consideration in the determination of 
atmospheric refraction. The following conditions 
specify the criteria for rejection of data. 

(1) Occurrence of dynamic solar activity. 

(2) R ado me not rotating. 

(3) Vertical axis maladjustment. 

(4) Periods of known abnormally large wind and 
temperature effects resulting from removal or 
addition of the radome. 

(5) Improper operation of certain components. 
The rejected interval was usually prior to complete 
failure or necessary readjustment of the component. 

(6) Transient tracking conditions; during acquisi- 
tion, reacquisition (after momentary power interrup- 
tion), radiometer operation, or tracking through an 
obstruction to the Line-of-sight. 

(7) Outer-loop offset-correction information inade- 
quate to compensate properly for this error source. 

(8) Data known to be bad at the time of observa- 
tion, and so labeled prior to computation of tracking 
error. 

The above data rejections eliminated observations 
which deviated significantly from the normal mode 
of operation, and thus tended to bring the observations 
into statistical control. 

There were approximately 40,000 observations 
made over a five-month interval which were selected 
and corrected for the three instrument errors and 
which were suitable for the analysis of atomospheric 
refraction. 

The observations of apparent altitude angle were 
made at 1-min intervals and hence were, in general, 
for nonintegral values of the observed angle. For 
this study, refraction was to be determined for the 
following observed altitude angles: 2, 4, 6, 8, 10, 14, 
18, 22, 26, 30, 35, 40, 45, 55, 60, and 65 deg. Thus, 
it was necessary to select the desired values of 
measured refraction from the entire collection of 
corrected observations. This was done by scrutiniz- 
ing the the data and listing all those observations 
whose observed altitude angle was within ±0.09 deg 
of the specified integral values. In cases where there 
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were more than one observation in this specified 
small interval, the one nearest the desired altitude 
was used. Surface refractivity, time of day, and 
status of the radome were listed along with measured 
refraction. 

Before proceeding with the averaging process at 
each desired altitude angle, it was necessary to de- 
termine whether or not there were differences in 
refraction at the same observed angle in the morning 
(a.m.) and the afternoon (p.m.). The data for 30 deg 
were examined, and all the a.m. — p.m. pairs with the 
radome on, were selected. Ten pairs were obtained 
and the difference between the a.m. sample average 
and the p.m. sample average was computed. This 
difference was tested for significance at the 5-percent 
level with the well-known ^-test for significance of 
means when the variances are known to be unequal. 
The test showed no significance, so it was repeated 
at 8 deg and 18 deg. The same result was obtained. 
It was felt that an inordinately large variance con- 
tributed to the lack of significance. Further exami- 
nation revealed that the pairs of samples obtained 
covered a time span of several months and hence 
there were large variations in surface refractivity. 
The test was repeated, using pairs taken only from 
the month of October. Angles of 10, 14, and 35 deg 
were added. The tests were now all significant ex- 
cept those at 8 deg which had only 4 pairs of observa- 
tions and at 14 deg which was very nearly significant. 
Thus, it was concluded that there was a real difference 
between the a.m. and p.m. values of refraction, and 
therefore, the averages must be obtained for both 
sets of altitude angles. 

In addition to the effect of a.m. and p.m., there 
was a question of whether or not the radome intro- 
duced a refractive effect of its own. Two different 
tests were employed to answer this question. The 
first was to simply test for a significant difference 
between a sample of observations with radome on 
and another sample of observations with radome off, 
when both samples were obtained for the p.m. alti- 
tude angle of 6 deg. In this case, the conclusion was 
that there is no significant difference in refraction 
with the radome on over that with the radome off. 
A second test considered a pair of observations (that 
is, radome on and off) as being the results of two 
measurements on the same test specimen when the 
specimen had been subjected to two different treat- 
ments; namely, radome on and radome off. It was 
attempted to make the individual specimens as 
homogeneous as possible by using p.m. conditions 
which had nearly (±1 Af-unit) the same surface 
refractivity. The differences between the refraction 
with radome off and radome on were determined for 
10 different specimens and the average difference was 
tested for statistical significance from zero. The 
result was that the average difference did not signfi- 
cantly differ from zero and therefore it was concluded 
that the radome did not introduce an important 
refractive effect of its own. 

These two preliminary tests provided the informa- 
tion necessary to decide how the observed refrac- 
tions should be grouped in order to obtain the best 



estimate of the average atmospheric refraction dur- 
ing the test period. That is, the observations must 
be separated into two groups depending on whether 
the desired altitude angle was obtained in the a.m. 
or in the p.m. It was not necessary to further 
separate the observations into radome on and radome 
off groups, hence they were lumped together. 

The final step in the computational process was to 
compute the mean, r, standard error of the mean, 
(Jt, and standard deviation, o>, of atmospheric re- 
fraction for each altitude angle, h . 

3. Discussion of Data and Results 

The results of computation of the statistics of 
refraction for each specified altitude angle are shown 
in table 1. 



Table 1. 



Statistics of measured atmospheric refraction for 
specified altitude angles at 1.85 cm 



deg 

6° a.m 

8° 

10° 

14° 

18° 

22° 

26° 

30° 

35° 

40° 

45° 

55° 

60° 

65° MP 

60°p.m 

55° 

45° 

40° 

35° 

30° 

26° 

22° 

18° 

14° 

10° 

8° 

6° 

4° 

2° 



T 


<Tt 


deg 


deg 


-0. 1702 


±0. 0048 


-. 1344 


±. 0042 


-. 1069 


±. 0027 


-. 0759 


±. 0016 


-.0608 


±. 0013 


-. 0480 


±.0011 


-.0403 


±. 0008 


-. 0348 


±. 0008 


-.0301 


±. 0007 


-. 0274 


±. 0009 


-. 0230 


±. 0012 


-.0171 


±. 0007 


-.0152 


±. 0009 


-.0112 


±. 0009 


-.0129 


±. 0007 


-.0153 


±. 0007 


-.0205 


±. 0007 


-.0237 


±. 0006 


-.0262 


±. 0006 


-.0328 


±. 0007 


-.0384 


±. 0008 


-.0465 


±. 0010 


-.0565 


±.0011 


-.0745 


±. 0014 


-. 1052 


±. 0019 


-. 1300 


±. 0021 


-. 1689 


±. 0030 


-.2342 


±. 0045 


-.3658 


±. 0073 



<Tt 


Number of 




observations 


deg 




11. 65X10-3 


6 


13.20X10-3 


10 


9. 11X10-3 


11 


6. 48X10-3 


17 


6. 18X10-3 


24 


5. 39X10-3 


25 


4.02X10-3 


24 


3. 72X10-3 


23 


3. 66X10-3 


25 


4. 25X10-3 


21 


5.01X10-3 


17 


2. 66X10-3 


13 


2.95X10-3 


10 


1. 90X10-3 


5 


2. 39X10-3 


11 


2. 70X10-3 


13 


2. 97X10-3 


21 


2. 92X10-3 


26 


3. 56X10-3 


33 


4. 15X10-3 


35 


5. 19X10-3 


39 


5. 10X10-3 


41 


6.40X10-3 


35 


8. 28X10-3 


34 


11. 30X10-3 


36 


13. 10X10-3 


39 


18.00X10-3 


36 


25. 80X10-3 


33 


40. 50X10-3 


31 



% 

2.8 
3.1 
2.5 
2.1 
2.1 
2.2 
1.9 
2.2 
2.3 
3.2 
5.2 
4.0 
5.9 
8.0 
5.4 
4.5 
3.4 
2.5 
2.2 
2.1 
2.0 
2.1 
1.9 
1.8 
1.8 
1.6 
1.7 
1.9 
1.9 



a =0-7/7X100. 

The table shows that these measurements of 
refraction are consistent with the form of the ex- 
pected variation of refraction. That is, at low 
angles, refraction is large (on the order of 0.5 deg) 
and it becomes very small at high-altitude angles, 
being minimum at meridian passage. Plots of the 
average refraction and standard error of estimate of 
the average versus observed altitude angle for a.m. 
and p.m. are shown in figures 1 and 2. 

The well-behaved smooth variation of r with alti- 
tude angle suggests that the observations were 
indeed measurements of refraction and that removal 
of extraneous effects has been rather successfully 
accomplished. Furthermore, it suggests that the 
radio sextant has been able to perform these measure- 
ments with high precision. 

The standard error of estimate of refraction re- 
mains just slightly less than 0.001 deg for all angles 
greater than 20 deg. Below 20 deg the standard 
error of estimate of the average increases with de- 
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creasing altitude angle. This undoubtedly results 
primarily from the increasing variability of refraction 
at low angles, and to a secondary extent, from the 
decreasing number of observations at very low angles. 

Table 1 and figures 3 and 4 show the standard 
deviation of the measured values of refraction for 
each altitude angle. The gross variation is similar 
to refraction in that the standard deviation is small 
at high-altitude angles and increases as altitude 
angle decreases. 

There are four factors which contribute to this 
increase of variability at low angles. The first 
involves the tracking system behavior with variation 
of loop gain. Tracking-servo random errors in- 
crease as the difference between the power received 
from the sun and the surrounding sky decreases. 
At low-altitude angles, the sun's radiation must pass 
through a greater distance of attenuating atmosphere 
than at high angles. This results in an appreciable 
reduction of power received from the sun at low 
angles. In addition, energy emitted from the 
atmosphere itself increases as the transmission path 
length increases. Although this second effect is 
very small, it may contribute to the less precise 
ability to distinguish between the sun and the sur- 
rounding sky at very low angles. The increasing 
path length, in effect, increases the tracking-servo 
random errors at low angles. 

The second factor is the rapidly changing meteor- 
ological conditions at the sunrise and sunset periods. 
At sunrise, insolation produces heating of the earth's 
surface and causes parcels of surface air to be heated 
and then to rise by convection. In addition to 
causing rapid variation in the surface refractivity, 
this heating from below also produces a convectively 
mixed lower atmosphere. As insolation decreases 
in the afternoon, convective mixing ceases, the sur- 
face temperature begins to decrease, and the atmos- 
phere begins to settle and stratify into more uni- 
formly varying layers. Near sunset, surface refrac- 
tivity is again changing appreciably, and the 
atmosphere is less homogeneous. 

The third factor contributing to increased varia- 
bility at low angles is the fact that the refractivity 
profile becomes more important at low angles. 
Thus, the state of the atmosphere above the surface 
contributes to the total refraction. Furthermore, 
as the altitude angle decreases, the ray path extends 
over greater distances from the observing station. 
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Figure 1. Average atmospheric refraction and standard 
of estimate obtained during a.m. 
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Figure 2. Average atmospheric refraction and standard 
of estimate obtained during p.m. 
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Figure 3. Standard deviation of refraction for a.m. observations. 
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Figure 4. Standard deviation of refraction for p.m. observations. 
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Thus, nonlocal variations in meteorology have an 
increasing effect on total refraction. 

The final contribution to variability of total re- 
fraction results from the inability to obtain meas- 
urements of refraction at exact integral values of 
observed altitude angle. Since the various values of 
observed angle were symmetrically distributed about 
the desired integral value, no systematic bias was 
introduced, but some of the variability in <j t results 
from this variability in h . The fact that h was 
truncated at ±0.09 deg from the desired value per- 
mits some estimate of the rms variability of the 
error in h , a A h , to be calculated in the usual man- 
ner by averaging the sum of squares of error. Thus 



°A» n — 






(h' —ho) 2 dho 



4! 

= 3 



VAh n 



0.09 



V3 1.73 



=0.052 deg, 



where h' Q is the desired integral value of observed 
altitude and A is the one-way maximum error (0.09 
deg). Using the slope of the r versus h curve, the 
variability in h can be converted into variability 
in r. Thus, 

The magnitude of the rate of change, {dr/dhol, can 
be found from the graph of refraction, which gives 
the maximum rate magnitude as 0.066 deg per degree 
at 2 deg altitude angle. Thus, the maximum esti- 
mated variability in refraction due to variability in 
selecting h is a' T = 0.0034 deg. This represents 8.4 
percent of the total variability in r at 2 deg altitude. 

In figure 5 are plotted the various values of the 
percent error of estimate of total refraction. The 
graph shows minimum percent error at the low alti- 
tude angles with a general increase to the maximum 
of 8 percent at 65 deg. This behavior arises from 
the fact that the standard error of estimate reaches 
a minimum value at high angles, but the total refrac- 
tion continues to decrease monotonically. 

The statistical tests described earlier, indicated 
that a.m. refraction was greater than p.m. refraction. 
Comparison of the two values for each altitude angle 
shown in table 1 verifies that this is consistently 
true. The differences are plotted as a function of 
altitude angle in figure 6. The dashed line indicates 
the average difference of 0.0025 deg. While this 
investigation did not attempt to determine why this 
is so, it is evidently true, and the following hypothesis 
is offered without verification. As a result of 
asymmetry in the diurnal temperature cycle, morn- 
ing hours are generally cooler than afternoon. Since 
refractivity is inversely proportional to temperature, 
the surface refractivity would be expected to be 
somewhat greater during the morning hours. Figure 
7 shows the average surface refractivity over a 



month for each hour of the day. The diurnal 
variation is clearly evident and shows that, on the 
average, a.m. surface refractivity in fact is greater 
than p.m. The greater surface refractivity and the 
more homogeneous mixing would tend to produce 
more refractive bending in the first few kilometers 
above the surface during the morning hours. 

Comparisons of measured radio refraction with 
standard optical refraction from the Nautical 
Almanac, 8.7-mm measured refraction, and that 
calculated by ray tracing from an assumed model 
atmosphere are shown in figure 8. This figure shows 
that the refraction measured by the radio sextant 
at 1.85 cm is consistently greater than any of the 
other determinations of atmospheric refraction. 
This is probably not a deficiency in the measure- 
ments described in this paper, but rather, is the 
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Figure 5. Percent error of estimate of total refraction. 
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result to be expected in view of the facts that the 
average surface refractivity during the sampling 
period was relatively high, and that 1.85-cm radio 
refractivity is higher than thai at optical wave- 
lengths. The latter results from the decreased 
dielectric polarization of the water vapor molecule. 

The difference between radio and standard optical 
refraction is seen to be appreciable even at high- 
altitude angles. Since the optical refraction was 
determined for standard conditions of temperature 
and pressure, it is believed that some reduction of 
this difference would be obtained if optical refrac- 
tivity were determined for the nonstandard con- 
ditions obtained during observation. Optical refrac- 
tivity does not include the contribution of water 
vapor, so a fundamental difference will exist no 
matter how well the optical refraction may be 
determined. Recent publications of the Nautical 
Almanac do have tables winch provide further 
correction to the standard optical refraction for 
nonstandard conditions. 

Calculated refraction gives the closest comparison 
with the measured results. The calculations were 
performed for a nonreal atmosphere with a surface 
refractivity of 318 TV-units, which is about 4.3 per- 
cent less than the average measured surface refrac- 
tivity during the experimental period. In other 
words, if the assumed model profile had a surface 
refractivity of 332 JV-units, the measured values of 
total refraction would then compare more favorably, 
especially at the 4 higher altitude angles. 

The second general feature of figure S is increasing 
discrepancy in all cases at angles below approxi- 
mately 20 degrees. This behavior is primarily a 
consequence of the difference in surface refractivity 
between the measured data and the optical and 
assumed models and, to a lesser extent, of the 
dependence, at low angles, of total refraction on the 
refractivity profile. That is, the total discrepancy 
at low angles is not accounted for by the difference 
in surface refractivity alone. A small portion of the 
low-angle discrepancy is attributed to difference in 
shape of the refractivity profiles for the measured 
data and the assumed models. 

Peculiar behavior of the 8.7-mm difference plot is 
believed due primarily to the assumption by Marner 
and Eingoen [1956] that change in refraction is 
negligible above about 30 deg. The inflection around 
10 deg probably results from systematically readjust- 
ing the tracking servo bias at 10 deg altitude angle. 

Theoretical considerations indicate that total re- 
fraction is a linear function of surface refractivity 
for reasonably high altitude angles [Bean and 
Cahoon, 1957]. To test the hypothesis that this 
relation holds for low-altitude angles, the typical 
results shown in figure 9 were obtained. At 8 deg 
altitude, the linear correlation coefficient was 0.87, 
which indicates a high degree of confidence in the 
linear hypothesis. Furthermore, the standard error 
of prediction is reduced from 0.014 to 0.007 by use 
of the linear predictor. There is little apparent 
curvature in the distribution of these sample points, 
which gives further support to the hypothesis. 
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Figure 8. Comparison of refraction differences (1.85 cm 
measured value minus specified value). 
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4. Conclusions 

A precise radiometric method of determining the 
average total atmospheric refraction at a 1.85-cm 
wavelength has been described. The maximum 
standard error of estimate of refraction was 0.0073 
deg at 2 deg altitude angle, and the minimum was 
0.0006 deg at 35 deg altitude angle. Expressed in 
terms of percentage of the mean, the error of esti- 
mate of the mean ranged from 8.0 percent to 1.6 
percent. Less than 8.4 percent of the total errors 
of measurement were due to the inability to make 
the observations precisely at the desired altitude 
angle. Errors due to this cause were maximum at 
the lowest altitude angles and decreased rapidly 
with increasing altitude angle. 

The specific determinations of refraction obtained 
during this investigation are valid only for August 
through December 1959 at Cedar Rapids, Iowa. 
Any extrapolation for other time or location must 
be used with care. If the average surface refrac- 
tivity for some other location and time is known, 
these results, could be used by modifying them in 
the ratio of that average value to 332 iV-units. 
Such a procedure should give a valid first approxi- 
mation to the average radio refraction for a new 
time and location. 

Probably the most significant conclusion which 
can be drawn from this investigation is the support 
which it gives to the various theoretical considera- 
tions which have been advanced. Notably among 
them is the analytical solution given by Smart [1956] 
and the linear relationship between surface refrac- 
tivity and total atmospheric refraction suggested by 
Bean and Cahoon [1957]. The success of these 
measurements also contributes considerable weight 
to the efficacy of the Smith and Weintraub [1953] 
equation for radio refractivity. 

This investigation has brought into somewhat 
sharper focus two particular results. They are the 
difference between a.m. and p.m. refraction at the 



same observed altitude angle, and the dispersion of 
radio refractivity. Both of these phenomena have 
theoretical bases and the present results suggest 
that further investigation would be warranted. 

The results reported here are strictly long-term 
average values and must not be interpreted as 
having any particular relation to the "instantaneous" 
refraction under other observing conditions. The 
results do, however, indicate the gross variational 
pattern and supply information on the range of 
variation to be expected throughout the fall months 
in the mid western United States . 
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